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We report the observation of the X(3823) in the process e + e“ —> 7r + 7r _ X'(3823) —> 7r + 7r _ 7Xci with a 
statistical significance of 6.2cr, in data samples at center-of-mass energies =4.230, 4.260, 4.360, 4.420 and 
4.600 GeV collected with the BESIII detector at the BEPCII electron positron collider. The measured mass of 
the X(3823) is (3821.7 ± 1.3 ± 0.7) MeV/c 2 , where the first error is statistical and the second systematic, 
and the width is less than 16 MeV at the 90% confidence level. The products of the Bom cross sections 
for e + e~ -4 7r + 7r”X(3823) and the branching ratio B[A'(3823) -4 7Xd. |C 2 ] are also measured. These 
measurements are in good agreement with the assignment of the X (3823) as the 'ffC’Df) charmonium state. 

PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Pq 


Since its discovery, charmonium - meson particles which 
contain a charm and an anti-charm quark - has been an excel¬ 
lent tool for probing Quantum Chromodynamics (QCD), the 
fundamental theory that describes the strong interactions be¬ 
tween quarks and gluons, in the non-perturbative (low-energy, 
long-distance effects) regime, and remains of high interest 
both experimentally and theoretically. All of the charmo¬ 
nium states with masses that are below the open-charm thresh¬ 
old have been firmly established 01 i; open-charm refers 
to mesons containing a charm quark (antiquark) and either 
an up or down antiquark (quark), such as D or 1). How¬ 
ever, the observation of the spectrum that are above the open- 
charm threshold remains unsettled. During the past decade, 
many new charmoniumlike states were discovered, such as 
the X(3872) H, the V(4260) 0 H and the Z c { 3900) 0- 
[7]. These states provide strong evidence for the existence of 
exotic hadron states 0], Although charged charmoniumlike 
states like the Z c (3900) provide convincing evidence for the 
existence of multi-quark states 0], it is more difficult to distin¬ 
guish neutral candidate exotic states from conventional char¬ 
monium. Moreover, the study of transitions between charmo- 
nium(like) states, such as the F(4260) -4- 7^(3872) fToh . is 
an important approach to probe their nature, and the connec¬ 
tions between them. Thus, a more complete understanding of 
the charmonium(like) spectroscopy and their relations is nec¬ 
essary and timely. 

The lightest charmonium state above the DD threshold is 
the V13770) Jl, which is currently identified as the 1 3 Di 
state [ly, the J = 1 member of the D -wave spin-triplet char¬ 
monium states. Until now there have been no definitive obser¬ 
vations of its two D -wave spin-triplet partner states, i.e., the 
1 3 D 2 and 1 3 D 3 . Phenomenological models predict that the 
1 3 D 2 charmonium state has large decay widths to yXd and 
7Xc2 0 . In 1994, the E705 experiment reported a candi¬ 
date for the 1 3 D 2 state with a mass of 3836 ± 13 MeV/c 2 
and a statistical significance of 2.8cr 0. Recently, the 
Belle Collaboration reported evidence for a narrow resonance 
X(3823) -4 7Xd in B meson decays with 3.8 ct significance 
and mass 3823.1 ± 1.8(stat) ± 0.7(syst) MeV/c 2 , and sug¬ 
gested that this is a good candidate for the 1 3 D 2 charmonium 


state 0. In the following, we denote the \ 3 D 2 state as ft 2 
and the ^>(3686) ['0(2S')] state as ft'. 

In this Letter, we report a search for the production of the if 2 
state via the process e + e _ -4 7r + 7r“X, using 4.67 fb -1 data 
collected with the BESIII detector operating at the BEPCII 
storage ring 0 at center-of-mass (CM) energies that range 
from yfs = 4.19 to 4.60 GeV lfl5tl . The ip 2 candidates 
are reconstructed in their 7 Xci and 7 Xc 2 decay modes, with 
Xd,c 2 —> 7 J/ip and J/f) -4 (l = e or /r). A GEANT4- 
based 0 Monte Carlo (MC) simulation software package is 
used to optimize event selection criteria, determine the detec¬ 
tion efficiency, and estimate the backgrounds. For the sig¬ 
nal process, we generate 40,000 e + e _ -4 7r + 7r _ X(3823) 
events at each CM energy indicated above, using an EVT- 
GEN lfl7| ] phase space model, with X(3823) -4 7 X^ 1 , C 2 - Ini¬ 
tial state radiation (ISR) is simulated with KKMC II1811 . where 
the Born cross section of e + e~ -4 7r + 7r _ 2T(3823) between 
4.1 and 4.6 GeV is assumed to follow the e + e - -4 Tt + n~ip' 
lineshape 0 - The maximum ISR photon energy is set to 
correspond to the 4.1 GeV/c 2 production threshold of the 
7 t + 7t - X( 3823) system. Final-State-Radiation is handled with 
PHOTOS |20|. 

Events with four charged tracks with zero net charge are se¬ 
lected as described in Ref. Q]. Showers identified as photon 
candidates must satisfy fiducial and shower quality as well as 
timing requirements as described in Ref. [21]. At least two 
good photon candidates in each event are required. To im¬ 
prove the momentum and energy resolution and to reduce the 
background, the event is subjected to a four-constraint (4C) 
kinematic fit to the hypothesis e + e _ -4 , that 

constrains the total four-momentum of the detected particles 
to the initial four-momentum of the colliding beams. The x 2 
of the kinematic fit is required to be less than 80 (with an ef¬ 
ficiency of about 95% for signal events). For multi-photon 
events, the two photons returning the smallest x 2 from the 4C 
fit are assigned to be the radiative photons. 

To reject radiative Bhabha and radiative dimuon 
( 7 e + e - / 7 /r + /r - ) backgrounds associated with photon 
conversion, the cosine of the opening angle of the pion-pair 
candidates is required to be less than 0.98. This restric- 
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tion removes almost all Bhabha and dimuon background 
events, with an efficiency loss that is less than 1% for 
signal events. The background from e + e“ —> r]J/ip with 
77 —> 7r + 7r _ 7r°/77r + 7r _ is effectively rejected by the in¬ 
variant mass requirement M('y'yTT + 'K~) > 0.57 GeV/c 2 . 
MC simulation shows that this requirement removes less 
than 1% of the signal events. In order to remove possible 
backgrounds from e + e _ — > 7 isr' 0 / ► 7isr J/ip, 

accompanied with a fake photon or a second ISR pho¬ 
ton, e + e~ r/ip' with 77 — > 77, and e + e _ —»■ 777//, 

the invariant mass of 7r + 7r ~J/ip is required to satisfy 
J/ip) — m(ip')\ > 6 MeV/c 2 ll22h . The sig¬ 
nal efficiency for the ip' mass window veto is 85% at 
i/s = 4.420 GeV and > 99% at other energies. 

After imposing the above requirements, there are clear 
J/ip peaks in the M(£ + £~) invariant mass distributions for 
the data. The J/ip mass window is defined as 3.08 < 
M(l + £~) < 3.13 GeV/c 2 . The mass resolution is de¬ 
termined to be 9 MeV/c 2 by MC simulation. In order to 
evaluate non -J/ip backgrounds, we define J/ip mass side¬ 
bands as 3.01 < M{£+£~) < 3.06 GeV/c 2 or 3.15 < 
M(£ + £~) < 3.20 GeV/c 2 , which are twice as wide as the 
signal region. The combination of the higher energy photon 
( 7 h) with the J/ip candidate is used to reconstruct Xci,c 2 sig¬ 
nals, while the lower one is assumed to originate from the 
X (3823) decay. We define the invariant mass range 3.490 < 
M(~/hJ/ ip) < 3.530 GeV/c 2 as the Xd signal region, and 
3.536 < M ( 7 //J/ 7 /)) < 3.576 GeV/c 2 as the Xc 2 signal re¬ 
gion [M(jbJ/i/j) = M("/h£ + £~) — M{£ + £~) + m(J/ip )]. 

To investigate the possible existence of resonances that may 
decay to 7Xci,c2, we examine two-dimensional scatter plots of 
M recoi i(n + TT-) versus M(j H J/ip). Here, M lecoi i(n + tt~ ) = 
y/(P e + e - — P-k+ ~ Pn-) 2 is the recoil mass of the 7 r + 7 r _ 
pair, where P e + e - and /V 1 are the 4-momenta of the ini¬ 
tial e + e _ system and the respectively. For this, we use 
the 7 r + 7 r~ momenta before the 4C fit correction because of 
the good resolution for low momentum pion tracks, as ob¬ 
served from MC simulation. Figure Q] shows M re coii( 7 r + 7 r _ ) 
versus M( 7 #J/t/>) for data at different energies, where 
e + e _ — y ir + n~ip' — > 7r + 7r _ 7X'ci,c2 signals are evident in 
almost all data sets. In addition, event accumulations near 
-^recoil( 7r+7I ' _ ) — 3.82 GeV/c 2 are evident in the Xci signal 
regions of the yfs = 4.36 and 4.42 GeV data sets. A scatter 
plot of all the data sets combined is shown in Fig. |l](f), where 
there is a distinct cluster of events near 3.82 GeV/c 2 (denoted 
hereafter as the AC(3823)) in the x c i signal region. 

The remaining backgrounds mainly come from e + e _ —y 
(rj'/'yuj)J/ip, with ( rj'/u)) —> 777 r + 7 r _ / 77 r + 7 r _ , and 
7 r + 7 r _ 7 r + 7 r _ ( 7 r°/ 77 ). The e + e _ (r)'/ / yui)J/ip back¬ 
grounds can be measured and simulated using the same data 
sets. The e + e _ —> n + n~n + n~ (tt°/' yy) mode can be evalu¬ 
ated with the J/ip mass sideband data. All these backgrounds 
are found to be small, and they produce flat contributions to 
the M reco ;i( 7 r + 7 r - ) mass distribution. There also might be 
e + e _ — > n + TT~ip' events with ip' —> r/J/ip and TT°n°J/ip, 
but such kind of events would not affect the ip' mass in the 
•Wrecoii( 7r+7r_ ) distribution. 



FIG. 1. Scatter plots of M Teco i\(n + iv ) vs. M(xhJ /ip) at (a) 
Vs =4.230, (b) 4.260, (c) 4.360, (d) 4.420, and (e) 4.600 GeV. The 
sum of all the data sets is shown in (f). In each plot, the vertical 
dashed red lines represent Xci (left two lines) and Xc 2 (right two 
lines) signal regions, and the horizontal lines represent the ip' mass 
range (bottom two lines) and 3.82 GeV (top line), respectively. 


An unbinned maximum likelihood fit to the M reco ii(7r + 7r _ ) 
invariant mass distribution is performed to extract the 
AC(3823) signal parameters. The signal shapes are repre¬ 
sented by MC-simulated ip' and X(3823) (with input mass of 
3.823 GeV/c 2 and a zero width) histograms, convolved with 
Gaussian functions with mean and width parameters left free 
in the fit to account for the mass and resolution difference 
between data and MC simulation, respectively. The back¬ 
ground is parameterized as a linear function, as indicated by 
the J/ip mass sideband data. The ip' signal is used to cal¬ 
ibrate the absolute mass scale and the resolution difference 
between data and simulation, which is expected to be similar 
for the V (3823) and ip'. A simultaneous fit with a common 
X (3823) mass is applied to the data sets with independent sig¬ 
nal yields at yfs = 4.230, 4.260, 4.360,4.420 and 4.600 GeV 
(data sets with small luminosities are merged to nearby data 
sets with larger luminosities), for the 7Xci an d 7Xc2 modes, 
respectively. 

Figure[2]shows the fit results, which return M[X(3823)] = 
M[X(3823)]i nput + yix ( 3823 ) ~ = 3821.7 ± 1.3 MeV/c 2 

for the 7Xci mode, where M[X(3823)]i npu t is the in¬ 
put X(3823) mass in MC simulation, Hx( 3823 ) = 1-9 ± 
1.3 MeV/c 2 and = 3.2 ± 0.6 MeV/c 2 are the mass shift 
values for X (3823) and ip' histograms from the fit. The 
fit yields 19 ± 5 X(3823) signal events in the 7Xci mode. 
The statistical significance of the X (3823) signal in the 7Xci 
mode is estimated to be 6.2cr by comparing the difference 
between the log-likelihood value (A(ln£) = 27.5) with or 
without X (3823) signal in the fit, and taking the change of 
the number of degrees of freedom (Andf = 6) into account. 
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FIG. 2 . Simultaneous fit to the M reco ii(7r + 7r _ ) distribution of 7X" c i 
events (left) and 7 x 72 events (right), respectively. Dots with error 
bars are data, red solid curves are total fit, dashed blue curves are 
background, and the green shaded histograms are J/0 mass side¬ 
band events. 
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FIG. 3. (a) The X(3823) scattering angle distribution for X(3823) 
signal events, the inset shows the corresponding M(n + n~) in¬ 
variant mass distribution per 20 MeV/c 2 bin; and (b) fit to the 
energy-dependent cross section of <r s [e + e - — > 7r + 7r _ A'(3823)] • 
B(X( 3823) —> 7Xci) with the V(4360) (red solid curve) and the 
0(4415) (blue dashed curve) lineshapes. Dots with error bars are 
data. The red solid (blue dashed) histogram in (a) is MC simulation 
with D -wave (S'-wave). 


and its value is found to be larger than 5.9 ct with various 
systematic checks. For the 7 Xc 2 mode, we do not observe 
an X (3823) signal and provide an upper limit on its produc¬ 
tion rate (Table HJ). The limited statistics preclude a measure¬ 
ment of the intrinsic width of X(3823). From a fit using a 
Breit-Wigner function (with a width parameter that is allowed 
to float) convolved with Gaussian resolution, we determine 
r[X(3823)] < 16 MeV at the 90% confidence level (C.L.) 
(including systematic errors). 

The A03823) is a candidate for the 0 a charmonium state 
with J PC = 2 lfl3tl . In the e + e _ — > 7T + 7r tp 2 process, 
the 7r + 7r - system is very likely to be dominated by S'-wave. 
Thus, a D- wave between the 7r + 7r _ system and 02 is ex¬ 
pected, with an angular distribution of 1 + cos 2 9 for 02 in 
the e + e _ CM frame. Figure [3] (a) shows the angular distri¬ 
bution (cos 9) of X (3823) signal events selected by requiring 
3.82 < M reco ii(7r + 7r - ) < 3.83 GeV/c 2 . The inset shows 
the corresponding invariant mass distribution per 

20 MeV/c 2 bin. A Kolmogorov [23] test to the angular dis¬ 
tribution gives the Kolmogorov statistic D 44 obs = 0.217 for 
the D- wave hypothesis and £>f 4 obs = 0.182 for the S'-wave 
hypotheses. Due to limited statistics, both hypothesis can be 
accepted (D B 4 obs , £>f 4 obs < £> 14 , 0.1 = 0.314) at the 90% 
C.L. 


The product of the Born-order cross section and the 
branching ratio of X (3823) -7 7Xci,c2 is calculated 

-7- 7r+7r-X(3823)] • £>[X (3823) -4> 


using a B [e + e 


7Xcl,c2j = 


N° 


' i-n| 2 ‘ 


-, where 1Y°-K 2 is the 


ber of A03823) 7 Xci,c 2 signal events obtained from 
a fit to the M reco ii(7r + 7r _ ) distribution, £j nt is the inte¬ 
grated luminosity, e is the detection efficiency, B c i jC 2 is 
the branching fraction of Xci,c2 -A 7J /0 — > and 

(1 + 6) is the radiative correction factor, which depends on 
the lineshape of e + e~ —> 7t + 7t - X( 3823). Since we ob¬ 
serve large cross sections at y/s = 4.360 and 4.420 GeV, 
> 7r + 7r _ J V(3823) cross section fol- 
y over the full energy range 


+ c“ 


we assume the e + e 
lows that of e + e~ 


of interest and use the e + e _ —> lineshape from 

published results jTf)] as input in the calculation of the ef¬ 
ficiency and radiative correction factor. The vacuum polar- 
is calculated from QED with 0.5% un- 


ization factor | 1 _ n p 
certainty [24]. The results of these measurements for the 
data sets with large luminosities at y/s = 4.230, 4.260, 
4.360, 4.420 and 4.600 GeV are listed in Table Q] Since at 
each single energy data the A03823) signal is not very sig¬ 
nificant, upper limits for production cross sections at the 
90% C.L. based on the Bayesian method are given [system¬ 
atic effects are included by convolving the X (3823) sig¬ 
nal events yield (n yield ) dependent likelihood curves with 
a Gaussian with mean value zero and standard deviation 
n yield • er sys , where cr sys is the systematic uncertainty of the 
efficiencies]. The corresponding production ratio of TZ^> = 
o- B [e+e~->7r+7i-~X(3823)]-B[X(3823)-»7Xei] 

<t b [e+ e~ — -Blip'— >-7Xci] 

= 4.360 and 4.420 GeV. 


is also calculated at 


We fit the energy-dependent cross sections of e + e _ —> 
tt + tt~X (3823) with the Y (4360) shape or the ■0(4415) shape 
with their resonance parameters fixed to the PDG values [2|]. 
Figure [3] (b) shows the fit results, which give D^ hs = 0.151 
for the V(4360) hypothesis (HI) and £>^ 2 bs = 0.169 for 
the t)>(4415) hypothesis (H2), based on the Kolmogorov test. 
Thus, we accept both the Y (4360) and the -0(4415) hypothe¬ 
ses (-^obs* D f 2 obs < d 5,o.i = 0.509) at the 90% C.L. 


The systematic uncertainties in the X (3823) mass measure¬ 
ment include those from the absolute mass scale, resolution, 
the parameterization of the X (3823) signal, and the back¬ 
ground shape. Since we use the 0' signal to calibrate the 
fit, we conservatively take the uncertainty of 0.6 MeV/c 2 in 
the calibration procedure as the systematic uncertainty due to 
the mass scale. The resolution difference between the data 
and MC simulation is also estimated by the 0 ' signal. Vary¬ 
ing the resolution parameter by ± 1 a, the mass difference in 
the fit is 0.2 MeV/c 2 , which is taken as the systematic un¬ 
certainty from resolution. In the X (3823) mass fit, a MC- 
simulated histogram with the width of X (3823) set to zero is 
used to parameterize the signal shape. We replace this his¬ 
togram with a simulated A03823) resonance with a width of 
1.7 MeV [13] and repeat the fit; the change in the mass for 
this fit, 0.2 MeV/c 2 , is taken as the systematic uncertainty due 
to the signal parameterization. Likewise, changes measured 
with a background shape from MC-simulated (r)'/'yu>)J/ip 
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TABLE I. Number of observed events ( N ° bs ), integrated luminosities (£) detection efficiency (e) for the A(3823) —> 7X' c i mode, 
radiative correction factor (1 ± S), vacuum polarization factor ( jjryn? ), measured Born cross section a B (e + e~ -A 7r + 7r _ X(3823)) times 

,6i(X(3823) — > 7 Xci) (erf • Si) and fi2(X(3823) —> 7Xc2) (erf • B 2 ), and measured Bom cross section a B (e + e~ — > (erf,) at 

different energies. Other data sets with lower luminosity are not listed. The numbers in the brackets correspond to the upper limit measurements 
at the 90% C.L. The relative ratio a W + e -^rr+n x(3823)]B(X(3823)-»7Xci) j s a j so ca i cu i a ted. The first errors are statistical, and the 

V [e+e~ —>-7T+7r — i^7Xcl) 

second systematic. 


Vs (GeV) 

C (pb- 1 ) 

AT 

obs 


e 

1 + <5 

1 

|i-np 

a B - 

Si (pb) 


erf • 

S 2 (pb: 

1 erf/ (pb) 

TZ^i 

4.230 

1092 

0 7+ 1 - 4 
u - 1 - 0.7 

« 

3.8) 

0.168 

0.755 

1.056 

0.12 t° 0 2 A 

± 

0.02 

« 

0.64) 


- 

34.1 ±8.1 ±4.7 

- 

4.260 

826 

1 1+ 1 - 8 
± * ± -1.2 

« 

4.6) 

0.178 

0.751 

1.054 

0.23l“; 3 ® 

± 

0.04 

« 

0.98) 


- 

25.9 ±8.1 ±3.6 

- 

4.360 

540 

3.9l 2 ;? 

« 

8.2) 

0.196 

0.795 

1.051 

1 in +0 - 64 

1 - 1U — 0.47 

± 

0.15 

(< 

2.27) 

(< 

1.92) 

58.6 ± 14.2 ±8.1 

0-20f q!q 

4.420 

1074 

7 c+3.6 

« 

13.4) 

0.145 

0.967 

1.053 

1 OQ+0-59 
i.z,o_ 0-4 6 

± 

0.17 

« 

2.19) 

« 

0.54) 

33.4 ± 7.8 ± 4.6 

o.39tg:£ 

4.600 

567 

i-9—i:? 

« 

5.4) 

0.157 

1.075 

1.055 

0.47tg;£ 

± 

0.07 

« 

1.32) 


- 

10.4±^ ± 1.5 

- 


events or a second-order polynomial indicate a systematic un¬ 
certainty associated with the background shape of 0.2 MeV/c 2 
in mass. Assuming that all the sources are independent, the 
total systematic uncertainty is calculated by adding the indi¬ 
vidual uncertainties in quadrature, resulting in 0.7 MeV/c 2 for 
the A (3823) mass measurement. For the A (3823) width, we 
measure the upper limits with the above systematic checks, 
and report the most conservative one. 

The systematic uncertainties in the cross section measure¬ 
ment mainly come from efficiencies, signal parameterization, 
background shape, decay model, radiative correction, and lu¬ 
minosity measurement. The luminosity is measured using 
Bhabha events, with an uncertainty of 1.0%. The uncer¬ 
tainty in the tracking efficiency for high momenta leptons is 
1.0% per track. Pions have momenta that range from 0.1 to 
0.6 GeV/c, and the momentum-weighted uncertainty is 1.0% 
per track. In this analysis, the radiative transition photons 
have energies from 0.3 to 0.5 GeV. Studies with a sample of 
J/ip -A pn events show that the uncertainty in the reconstruc¬ 
tion efficiency for photons in this energy range is less than 
1.0%. 

The same sources of signal parameterization and back¬ 
ground shape as discussed in the systematic uncertainty of 
X (3823) mass measurement would contribute 4.0% and 8.8% 
differences in A(3823) signal events yields, which are taken 
as systematic uncertainties in the cross section measurement. 
Since the X (3823) is a candidate for the ip 2 charmonium state, 
we try to model the e + e _ — > ir + tt~X (3823) process with a 
D -wave in the MC simulation. The efficiency difference be¬ 
tween D -wave model and three-body phase space is 3.8%, 
which is quoted as the systematic uncertainty for the decay 
model. The e + e _ -A 7 t + 7t _ A( 3823) lineshape affects the 
radiative correction factor and detection efficiency. The radia¬ 
tor function is calculated from QED with 0.5% precision [25). 
As discussed above, both Y (4360) lineshapes fl9L l26ll and 
the ^(4415) lineshape describe the cross section of e + e _ —y 
7t + 7t - A( 3823) reasonably well. We take the difference for 
(1 + 5) ■ e between Y (4360) lineshapes and the ^(4415) line- 
shape as its systematic uncertainty, which is 6.5%. 

Since the event topology in this analysis is quite similar 


to e + e —> 77T + 7r J!ijj we use the same systematic 

uncertainties for the kinematic fit (1.5%) and the J/ip mass 
window (1.6%). The uncertainties on the branching ratios for 
Xci,c 2 —> 7 JM (3.6%) and J/tj) -A £ + £~ (0.6%) are taken 
from the PDG yj- The uncertainty fromMC statistics is 0.3%. 
The efficiencies for other selection criteria, the trigger simu¬ 
lation [27], the event-start-time determination, and the final- 
state-radiation simulation are very high (> 99%), and their 
systematic uncertainties are estimated to be less than 1%. 

Assuming that all the systematic uncertainty sources are in¬ 
dependent, we add all of them in quadrature. The total system¬ 
atic uncertainty in the cross section measurements is estimated 
to be 13.8%. 

In summary, we observe a narrow resonance, A (3823), 
through the process e + e _ -A 7r + 7r _ A(3823) with a statistical 
significance of 6.2er. The measured mass of the A(3823) is 
(3821.7±1.3±0.7) MeV/c 2 , where the first error is statistical 
and the second systematic, and the width is less than 16 MeV 
at the 90% C.L. Our measurement agrees well with the val¬ 
ues found by Belle [iTU. The production cross sections of 
a B (e+e~ -A 7r+7r-A(3823)) • B(X( 3823) -a 7Xd,7Xc 2 ) 
are also measured at y/s = 4.230, 4.260, 4.360, 4.420, and 
4.600 GeV. 

The A(3823) resonance is a good candidate for the 
i/j(1 3 D2) charmonium state. According to potential mod¬ 
els d, the D- wave charmonium states are expected to be 
within a mass range of 3.82 to 3.85 GeV. Among these, the 
1 1 D ‘2 -a 7Xci transition is forbidden due to C-parity con¬ 
servation, and the amplitude for 1 3 Z?3 -a 7Xci is expected 
to be small l28ll . The mass of ^(l 3 !^) is in the 3.810 ~ 
3.840 GeV/c 2 range that is expected for several phenomeno¬ 
logical calculations f29il . In this case, the mass of •?/’(1 3 Z7 2 ) 
is above the DD threshold but below the DD* threshold. 
Since ^(l 3 !^) -A DD violates parity, the t/>(l 3 £> 2 ) is 
expected to be narrow, in agreement with our observation, 
and ^(lflA) -A 7Xci is expected to be a dominant decay 
mode 1291130 1. From our cross section measurement, the ra¬ 
tio b[y( 3823 )^> 7 x C i] < (wh ere systematic uncertainties 
cancel) at the 90% C.L. is obtained, which also agrees with 
expectations for the ^(l 3 !^) state & 
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